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Abstract. 

A generalization of the Zernike circle polynomials for expansion of functions 
vanishing outside the unit disk is given. These generalized Zernike functions 
have the form = R^'^'^p) exp{M), < p < 1, < ^ < 2tt, 

and vanish for p > 1, where n and m are integers such that n — \m\ is non- 
negative and even. The radial parts are — as p t 1 which a 
is a real parameter > —1. The Z™'" are orthogonal on the unit disk with 
respect to the weight function (1 — p^)~°, 0<p<l. The Fourier transform 
of Zl^'°' can be expressed explicitly in terms of (generalized) Jinc functions 
J„+a+i(27rr)/(27rr)""'"-^ and exhibits a decay behaviour r~"~^/^ as r — )■ oo. 
This explicit result for the Fourier transform generalizes the basic identity in 
the classical Nijboer- Zernike theory, case a = 0, of optical diffraction. The 
generalized Zernike functions were considered by Tango in 1977 and a version 
of the result on the Fourier transform is presented, however, in somewhat dif- 
ferent form and without a detailed proof. The new functions accommodate 
the solution of various forward and inverse problems in diffraction theory 
and related fields in which functions vanishing outside the unit disk with 
prescribed behaviour at the edge of the disk are involved. Furthermore, the 
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explicit result for the Fourier transform of the allows formulation and 
solution of design problems, involving functions vanishing outside the unit 
disk whose Fourier transforms and their decay are prescribed, on the level 
of expansion coefficients with respect to the 2'^'". Cormack's result for the 
Radon transform of the classical circle polynomials admits an explicit gen- 
eralization in which the Gegenbauer polynomials C""^^ appear. This result 
can be used to express the radial parts as Fourier coefficients of 

C"^^ (p cos '&) and to devise a computation scheme of the DCT-type for fast 
and reliable computation of the i?™'"(p). 

In recent years, several (semi-) analytic computation schemes for forward 
and inverse problems in optical (ENZ) and acoustic (ANZ) diffraction theory 
have been developed for the classical case a — 0. Many of these schemes 
admit a generalization to the cases that a ^ 0. For instances that this is not 
the connection formula expressing the generalized Zernike functions 

as a linear combination of the classical ones can be used. Various acoustic 
quantities that can be expressed via King's integral (involving the Hankel 
transform of radially symmetric functions) can be expressed for the general 
case that a 7^ in a form that generalizes the form that holds for the case 
a = 0. Furthermore, an inverse problem of estimating a velocity profile 
on a circular, baffied piston can be solved in terms of expansion coefficients 
from near-field data using Weyl's representation result for spherical waves. 
Finally, the new functions, case a — ±1/2, are compared with certain non- 
orthogonal trial functions as used by Streng and by Mellow and Karkkainen 
with respect to their ability in solving certain design problems in acoustic 
radiation with boundary conditions on a disk. 
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1 Introduction 



In optical and acoustic radiation with boundary conditions given on a 
disk, the choice of basis functions to represent the boundary conditions is of 
great importance. Here one can consider the following requirements: 

A. the basis functions should be appropriate given the physical background, 

B. the basis functions should be effective and accurate in their ability of 
representing functions vanishing outside the disk and arising in the par- 
ticular physical context, 

C. the basis functions should have convenient analytic results for transfor- 
mations that arise in a natural way in the given physical context. 

For the case of optical diffraction with a circular pupil having both phase 
and amplitude non-uniformities, Zernike [1] introduced in 1934 his circle 
polynomials, denoted here as ^^(p, "&) with radial variable p > and angular 
variable vanishing for p > 1, which find nowadays wide-spread application 
in fields like optical engineering and lithography [2]-[7], astronomy [8j-|iIU] 
and ophthalmology |TT]-[T3]. The Zernike circle polynomials are given for 
integer n and m such that n — \m\ is even and non- negative as 

Z-(p,^) = /^^l(p)e^™^ 0<p<l, 0<^<27r, 

= , P>1, (1) 

where the radial polynomials /^IT' are given by 

i?H(p)=pHpfjH)(2/-l), (2) 



with -P^" the general Jacobi polynomial as given in [T3], Ch. 22, pL5], Ch. 4 
and [16], Ch. 5, §4. The circle polynomials were investigated by Bhatia and 
Wolf [17] with respect to their appropriateness for use in optical diffraction 
theory, see issue A above, and were shown to arise more or less uniquely 
as orthogonal functions with polynomial radial dependence satisfying form 
invariance under rotations of the unit disk. 
The orthogonality condition, 

1 27r 

J J Z„7(p,^)(Z-/(p,^)rprfprft9 = -^5^,^,5„,„, , (3) 
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where n is either one of rii and n2 at the right-hand side, together with 
the completeness, see [IS], App. VII, end of Sec. I, guarantees effective and 
accurate representation of square integrable functions on the unit disk in 
terms of their expansion coefficients with respect to the Z^, see issue B above. 
The Zernike circle polynomials, notably those of azimuthal order m = 0, were 
considered recently by Aarts and Janssen [TS]^[22] for use in solving a variety 
of forward and inverse problems in acoustic radiation. This raised in the 
acoustic community p3] the question how the circle polynomials with m = 
compare to other sets of non-polynomial, radially symmetric, orthogonal 
functions on the disk. It was shown in [23] that the expansion coefficients, 
when using circle polynomials, properly reflect smoothness of the functions 
to be expanded in terms of decay and that they compare favourably in this 
respect with the expansion coefficients that occur when using orthogonal 
Bessel series expansions; the latter expansions are sometimes used both in 
the acoustic and the optical domain, see [25] [27]. 

The present paper focuses on analytic properties, see issue C above, of 
basis functions, and in Sec. H] we present a number of such properties for 
the set of Zernike circle polynomials. In Sec. [3] a generalization of the set of 
Zernike circle polynomials is introduced, viz. the set of functions 

Zr"(p,t9) = (l-p2)"pHpKkn|)(2p2_i)g.rn. ^ 0<p<l, 

2 

= , p>l, (4) 

with p, ■}) and n, m as before, see ([T]), and parameter a > —1. These 
new Zernike-type functions arise sometimes more naturally in certain phys- 
ical problems and can be more convenient when solving inverse problems in 
diffraction theory since the decay in the Fourier domain can be controlled by 
choosing the parameter a appropriately. While issue A above should guide 
the choice of a, the matter of orthogonality and completeness is settled as 
in the classical case a = 0. A more involved question is what becomes of 
the analytic properties, noted for the classical circle polynomials in Sec. [2], 
when a 7^ 0. A major part of the present paper is concerned with answering 
this question, and this yields generalization of many of the results holding 
for the case a = and about which more specifics will be given at the end 
of Sec. [3] where the basic properties of the generalized Zernike functions are 
presented. 
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2 Analytic results for the classical circle poly- 
nomials 



The Zernike circle polynomials were used by Nijboer in his 1942 thesis 
for the computation of the point-spread functions in near best-focus planes 
pertaining to circular optical systems of low-to-medium numerical aperture. 
In that case, starting from a non-uniform pupil function 

P(p, ^) = A{p, ^) e'^^P'"^^ , < p < 1 , < ^9 < 27r , (5) 

the point-spread function U{r,ip ; f) at defocus plane / with polar coordi- 
nates X + iy = r e^'^ is given in accordance with well established practices in 
Fourier optics as 

1 27r 



f/(r,(^;/)= / / e'^''\^'''f"^'°'^^-'^^ P{p,^) pdpd^ . (6) 








It was discovered by Zernike [T] that the Fourier transform of the circle 
polynomials, 



1 2n 



(J-Zr)(r,^) = j j e2-^''^™^(''-^)zr(p,^9)pc^prf^ (7) 



has the closed-form result 

(J^Z:^){r,y,) = 27ie^-^^^^±^^e'"^^ , r>0, 0<y.<27r, (8) 

where Jn+i is the Bessel function of the first kind and order n + 1. Thus, one 
expands the pupil function P as 

P(p,^) = ^/3;r^r(P,^) , 0<p<l, 0<^9<27r, (9) 

n,m 

and uses the result in ([HD to compute the point-spread function f/ in (jH]) in 
terms of the expansion coefficients In the case of best focus, / = 0, this 
gives a computation result for the point-spread function immediately. For 
small values of /, say |/| < 1, one expands the focal factor, 

e'^^' = l + ^fp'-lfp'-... (10) 

and spends some additional effort to write functions p^'^ Z™(p, t?) as a k- 
terms linear combination of circle polynomials with upper index m to which 
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([7]) applies. Also see [H], Ch. 9, Sees. 2-4. These are the main features of 
the classical Nijboer-Zernike theory for computation of optical point-spread 
functions in the presence of aberrations. 

In [29] Janssen has computed the point-spread function U^{r, ip ; /) per- 
taining to a single term in the form 

f/-(r,¥,; /) = 27r^H^H(^j)e^-^ , (11) 

where 

1 

Vt^ {r,f) = I e'f^' i?!™! (p) J|„| (27rpr) p dp = 



Z=0 j=0 

with explicitly given (p = | — |^|))- This result has led to what is 
called the extended Nijboer-Zernike (ENZ) theory of forward and inverse 
computation for optical aberrations. The theory in its present form allows 
computation of point-spread functions for general / and for high-NA optical 
systems, including polarization and birefringence, as well as for multi-layer 
systems. See [SOI^ISI] and [3S] for an overview. Furthermore, it provides a 
framework for estimating pupil functions P, in terms of expansion coefficients 
/3, from measured data |f/p of the intensity point-spread function in the focal 
region. The latter inverse problem has the basic assumption that P deviates 
only mildly from being constant so that the term with m = n = in ([2]) 
dominates the totality of all other terms. The theoretical intensity point- 
spread function can then, with modest error, be linearized around the leading 
term |/3o t/gP; ^^i^ leads, via a matching procedure with the measured 
data in the focal region, to a first estimate of the unknown coefficients. This 
procedure is made iterative by incorporating the totality of all deleted small 
cross-terms (involving the (m, n) ^ (0, 0), quadratically) in the matching 
procedure using the estimate of the /3's from the previous steps. In practice 
one finds that pupil functions P deviating from being constant by as much 
as 2.5 times the diffraction limit can be retrieved. See [32], [33], [35], [36], 
[37] for more details. 

The result in ([H]) is one evidence of computational appropriateness of the 
circle polynomials for forward and inverse problems, but there are others. 
Many of these are based on the basic NZ-result in (|8]). In |38]-[39], Cormack 
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used this result to calculate the Radon transform 7^^, 

7e-(r,^)= I Z!:{u,f,)dl (13) 

of Z™, with integration along the line u cos ip + fJ^ sin = t in the plane with 
r > and ip £ [0, 27r) and where (i^, /x) = (pcos-f^jpsin-;?). The result is that 

^n(^,V^) = ^(l-r')'/'^/n(r)e^'"^ 0<r<l, < < 27r , (14) 
n + 1 

where f/„ is the Chebyshev polynomial of the second kind and degree n, |14] . 
Ch. 22. On this explicit form of the Radon transform of Z™, Cormack based 
a method for estimating a function on the disk from its Radon transform 
by estimating its Zernike expansion coefficients through matching. (In 1979, 
Cormack was awarded the Nobel prize in medicine, together with Hounsfield, 
for their work in computerized tomography.) Cormack's result was used by 
Dirksen and Janssen [30] to find the integral representation 

R^[p) = — f Unip COS COS m'&d^ (15) 

27T J 


(integer m, n > 0) that displays, for any p > 0, the value of the radial part at 
p in the form of the Fourier coefficient of a trigonometric polynomial. This 
formula ( !T5|) can be discretized, error free when more that m + n equidistant 
points in [0, 2n] are used, and this yields a scheme of the DCT-type for 
computation of i?™(p). 

A further consequence of the basic NZ-result (|8]) is the theory of shifted- 
and-scaled circle polynomials developed in [4^. For given a > 0, 6 > 
with a + b < 1, there are developed explicit expressions for the coefficients 
-^rm™ (a, b) in the expansion 

Z;^((a + 6p'cos?9',6p'sin?9')) = Yl ^^'^'(a, ^^'(p' cos?9, p' sin ?9') . (16) 

n' ,m' 

This generalizes the result in [12], 

Rn{ep)= Yl {Rt{e)-Ri+\s))K,{p) , (17) 

n'=r?i(2)n 

on the Zernike expansion of scaled circle polynomials (m > 0). 
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The basic NZ-result is also useful for the computation of various acoustic 
quantities that arise from a circular piston in a planar baffle. The complex 
amplitude p{r, k) of the sound pressure at the field point r = {x,y, z), z > 0, 
in front of the baffle plane z = due to a harmonic excitation exp{iu)t) with 
wave number k = cu/c, with c the speed of sound, is given by Rayleigh's 
integral and King's integral as 



oo 



^__^(„2_^2)l/2 



^Po'^^ I (^2 _ p)i/2 M^u) ^(^) ^ du . (18) 







Here po is the density of the medium, v{a) is a non-uniform velocity profile 
assumed to depend on the radial variable a = (z/^ + js^Y^"^ on the piston 
surface S with center and radius a, and r' is the distance from the field 
point r to the point (i/, yU, 0) on 5*. Furthermore, w = (x^ + y'^Y^'^ is the 
distance from the field point r to the 2;-axis, the root (m^ — k'^Y^'^ has the 
value i y/k^ — u"^ and y/u"^ ~ k^ for < u < k and u > k, respectively with 
non-negative in both cases, and V{u) is the Hankel transform, 

a 

V{u)= / Jo{u(r)v{(r)ada , u>0, (19) 



of order of v{a). 

In [TP], the on-axis pressure P2i{r = {0,0,z),k) due to v{a) = R^i^crja), 
< cr < a, was shown from Rayleigh's integral and a special result on 
spherical Bessel functions to be given as 

P2Kr = (0,0,^),A;) = \poc{ka)\-l)'ji{kr^)hf\kr+) , (20) 

where ji and hf"* are spherical Bessel functions, see [14j, Ch. 10 and in par- 
ticular 10.1.45-46, and r± = ^ [(r^ + a?Y^'^ ±r]. This result was used in [T^] 
for estimating a velocity profile from on-axis pressure data on the level of 
expansion coefficients with respect to radially symmetric circle polynomials. 
In [20], the King integral for the pressure is employed to express the 

pressure p((l,0,0),fc) at the edge, the reaction force J pdS and the total 

s 

radiated power J p(0) f* (a) (iS* in integral form. Expanding f (a) into radially 
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symmetric circle polynomials and using the basic NZ-result for an explicit 
expression of the Hankel transform V{u) in flT^ . this gives rise to integrals 



oo 



oo 



/ 



(m2 - A;2)l/2 



du , 



/ 



(u2 - A;2)i/2y 



du 



(21) 











with integer m, n > of same parity. These integrals have been evaluated as 
a power series in ka in [20] . 

In [21], the problem of sound radiation from a flexible spherical cap on a 
rigid sphere is considered. The scaling theory of Zernike circle polynomials, 
appropriately warped so as to account for the spherical geometry of the 
problem, is used to bring the standard solution of the Helmholtz equation 
with axially symmetric boundary data in a semi-analytic form per warped 
Zernike term. This gives rise to a coefficient-based solution of the inverse 
problem of estimating an axially symmetric velocity profile on the cap from 
measured pressure data in the space around the sphere. 

Finally, returning to baffled-piston radiation, in |22] the impulse response 
h{r, t), t > 0, at a field point r = (x, y, z) in front of the baffle is considered. 
This h{r,t) is obtained as a Fourier inversion integral with respect to wave 
number k of the velocity potential ip{r,k) = {ipock)~^ p{r, k), which can be 
evaluated, via Rayleigh's integral, as an integral of the velocity profile v over 
all points in the baffle plane at a common distance (c^t^ — z^)^/^ from the field 
point. For the case that v{a) = R2i{cr/a), < a < a, the latter integral can 
be evaluated explicitly using the addition theorem for Legendre polynomials. 
This explicit result can be used to compute impulse responses for general 
V by expanding such a v into radially symmetric circle functions. And, of 
course, an inverse problem, with a coefficient-based solution, can again be 
formulated. 

3 Generalized Zernike circle functions 

The Zernike circle polynomials, see ([T]), all have modulus 1 at the edge 
p = 1 of the unit disk while in both Optics and Acoustics it frequently occurs 
that the non-uniformity behaves differently towards the edge of the disk. 

In optical design, it is often desirable to have a pupil function P whose 
point-spread function U, see ([6]), decays relatively fast outside the focal region 
where the specifications of the designer are to be met. The Zernike circle 
polynomials are discontinuous at the edge p = 1, and the corresponding 
point-spread functions U^{r, ip ; /) have poor decay, like r~^/'^ as r — )• oo. 
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In the ENZ point-spread functions, see Sec. [21 the key assumption is that 
the pupil function's deviation from being constant is not large. When \P\ 
gets small at the edge, which happens, for instance, when the source is a 
pinhole with a positive diameter while the objective lens has a large NA, this 
basic assumption is not met. In such aberration retrieval with the 

ENZ method may become cumbersome, even in its iterative version. 

In the theory of acoustic radiation from a baffled, planar piston with 
radially symmetric boundary conditions, Streng [13] follows the approach of 
Bouwkamp [?^J in solving the Helmholtz equation for this case and postulates 
normalized pressure functions on the disk of the form 

oo 

a^l - p')"^'/\ 0<p<l. (22) 

n=0 

In the light of the later developments of this paper, it is interesting to note 
here that Bouwkamp himself prefers expansions that involve the functions 

(1 - p'y/' P/^/^'°)(2p2 - 1) = (-1)' P2m((l - P'V^') , < p < 1 . (23) 

Similarly, Mellow has velocity profiles on the disk of the form (normalized 
to the unit disk) 

oo 

v{p) = J2bn{l-p')''-'^\ 0<p<l. (24) 

n=0 

In |15] and [l6]-[17j, the postulates fl22l) . fl2^ are used to solve design prob- 
lems in acoustic radiation from a membrane in a circular disk in terms of the 
expansion coefficients a„, bn- In principle, these design problems could also 
be solved when radially symmetric circle polynomials instead of (1 — p2)"='=i/2 
were used as trial functions, but this would become cumbersome, the func- 
tions (1 — p^)^^/^ themselves already having a poorly convergent expansion 
with respect to the i?2;(p)- 

In the three instances just discussed, it would be quite helpful when the 
system of Zernike circle polynomials would be replaced by systems whose 
members exhibit an appropriate behaviour at the edge of the disk. These 
new systems should satisfy the requirements A, B and C mentioned in the 
beginning of Sec. [TJ In this paper, the set of functions 

Z™'"(p, ^) = (1 - p'r P'"' PtS\2p' - 1) e^'"" , 0<p<l,0<^9<27r 



= , P > 1 , (25) 
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is considered for this purpose. In f l2^ . the parameter a > — 1, and n and m 

Die 

k 



are integers such that n — \m\ is even and non-negative, and the Pjf''^'^ are 



Jacobi polynomials as before. Observe that the radial parts 

i?lrl'"(p) = (1 - p2)>l'"l PtS^^P' - 1) ' < p < 1 , (26) 

2 

depend non-polynomially on p, unless a = 0, 1, ... . The radial parts occur 
essentially in Tango [IB], Sec. 1, but there the factor (1 — p^)" is replaced by 
(1 — p^)°/^ with somewhat different restrictions on a and |m|. There is also 
a relation with the disk polynomials as they occur, for instance, in the work 
of Koornwinder [l9] : 

Zr"(p,^) = (1 ^{pe'") , (27) 

2 ' 2 

with 

DUpe'") = P"^-" - 1) e^^'^-'^" (28) 

for k,l = 0,1,.... Thus, the disk polynomials omit the factor (1 — p^)" 
altogether. 

Further definitions used in this paper are 

n — \m\ n + Iml , , 

P=^. -J^^ (29) 

for integers n and m such that n — \m\ is even and non- negative, and the 
generalized Pochhammer symbol 



By orthogonality of the Plf'^\x) with respect to the weight function 
'1 - a;)°(l + x)^ on [-1, 1], there holds 



1 27r 

1 

TT 




l-p2)-"Z--"(p,t9)(Z-/'°(p,^9))*pdpd^ 







(p ~l~ l)a ^•m\m2 ^rain2 



(31) 



(p + |m| + l)a ni + a + 1 
We note that for any integer = 0, 1, ... 

p\Mpi^){2p' -l)e''^^ (32) 

2 
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with n = 0, 1, and m = —n, —n + 2, n are | (A^ + 1)(A^ + 2) linearly 
independent functions of the form 'Yliij '^ij^V'' = pcosd, fi = psint?) 
where the summation is over integer i,j>0 with i + j < N. Therefore, by 
Weierstrass theorem the functions in ( 132|) are complete. 

This paper focuses on establishing versions for the generalized Zernike 
functions of the analytic results that were presented for the classical circle 
polynomials in Sec. HI Thus, in Sec. SJ the Fourier transform of is 
computed in terms of Bessel functions and a Weber-Schafheitlin representa- 
tion of the radial functions -R™'" is given. In Sec. O the Radon transform 
of Z™'" is computed in terms of the Gegenbauer polynomials C""^^, and a 
representation of the radial functions as Fourier coefficients of the periodic 
function C""^^(p cos < < 27r, for a fixed p G [0,1) is given. Thus a 
computation scheme of the DCT-type for the radial functions arises. Then, 
in Sec. El a scaling result is given. This result is somewhat more awkward 
than in the classical case since the radial functions i?^'°(p) have restrictions 
on their behaviour as p f 1 while the scaled radial functions R^'"{ep), to 
be considered with < e < 1, do not. Next, in Sec. [TJ the new Zernike 
functions are expanded in terms of the classical circle polynomials, with an 
explicit expression for the expansion coefficients. This makes it possible to 
transfer all forward computation results from the ENZ theory and from the 
acoustic Nijboer- Zernike (ANZ) theory for the classical case to the new set- 
ting in a semi-analytic form. This is useful for those cases that a closed 
form or a simple semi-analytic form is not available or awkward to find in 
the new setting. In Sees. [Hl-fTOl the focus is on generalizing the results from 
the ANZ theory, as developed in ^19j-^22j to the more general setting. This 
gives rise in Sec. [8] to a power series representation of the basic integrals that 
occur when various acoustic quantities are computed from King's integral 
for the sound pressure (in the case of baffied-piston radiation). An inverse 
problem, in which the velocity profile is estimated in terms of its expansion 
coefficients from near-field measurements via Weyl's formula, is considered 
in Sec. [91 Finally, in Sec. [TOl the trial functions as used by Streng and Mel- 
low, following Bouwkamp's solution of the diffraction problem for a circular 
aperture, are compared with the Z'^"~^^^'^ . In [4^, Mellow and Karkkainen 
consider radiation from a disk with concentric rings, and for this an indefinite 
integral involving the radial functions I^2t^^'^{p) is required. This integral is 
computed in closed form. 
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4 Fourier transform of generalized Zernike cir- 
cle functions 



In this section, the 2D Fourier transform of is computed. It is 

convenient to write here 

Z-'"(z/,/i)^Z-'"(p,.9) , (33) 
where u + ifi = pe^^ with p > and < < 2tt. 

Theorem 4.1. For a > — 1 and integer n, m such that p = — \m\) 
is a non-negative integer, there holds 



= 2n 2"(p + 1). ^""[^^l^ll^ , (34) 

where x + iy = r e'^'^ with r > and < < 27c. Furthermore, 
1 

I i?ri'"(p) JH(2vrpr)prfp= (-l)^2"(p + l).:^^g^ , (35) 



and 

oo 

i?l™l'"(p) = (-l)P 2°(p + 1), j Jn+a+l{t)J\m\{pt) Q < p < 1 . (36) 



Proof. First consider the case that m > 0. It follows from [13], 11.4.33 
(Weber-Schafheitlin integral), with 

p = n+p+1 = m + 2p+a + l , a = l, z/ = m, 6 = pG [0,1), \ = a (37) 
that 

oo 

Jn+a+l(t) Jm{pt) 



■dt 







p™r(m + p + l) ^, . . 2^ 

{—p — a,m + p + 1 ; m + 1; p J 



2"r(m + l)r(p + a + 1) 

(38) 
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Next from [H], 15.3.3 with 

a = m + l+ p + a, b = —p , c = m + l, z = (39) 
it follows that 

F{—p — a, m + p + 1 ; m + 1 ] p^) = 

= {1 - p'^y F{-p,m + l + a + p; m + 1; p'^) = 

= (1 - p'r T-XTV ^^"^(1 - V) , (40) 

where in the last step [H], 15.4.6 with n = p, a = m, I3 = a and z = p^ has 
been used. Therefore, for general integer m, |m| < n, 

Jn+a+l{t) J\m\{pt) _ 



P 



(1 -/)"p(°'H)(2p^ - 1) >^ /gH'"(p) 



2-(p + l)," V P / 2"(p+l), 

(41) 

where pI'''^\-x) = {-if P^^'"\-x) and the definition ^ of i?l"''"(p) 
have been used. This estabhshes (l36|) . 
Next, 

oo 27r 

(27rr)"+i ^ 








oo 



(-l)IH /• JHW^n+«+l(t)^,„,^, 



ini'd 



27r / 



cite*"^'^ , (42) 



where, subsequently, use is made of 

vx + py = pr cos(i9 — ip) , (43) 
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dd = J„(-t) = (-Ol'"! J|„|(t) , (44) 



see [13], Sec. 9.1, and where the substitution t = 27rr has been used in the last 
step in fH2]) . Therefore, using fHTl) . the definitions fl2^ - fl2Bl) and = |m|+2p, 
it follows that 

(27rr)"+i ^ 

'r"(p,^). (45) 



2n 2"(p+l)« 

Then flM|l follows by 2D Fourier inversion. Now also fl35l) follows using the 
definitions (125D-(I26D in §^ and proceeding as in fH^jl - fHij) . 

Notes. 

1. The result in ( l34l) generalizes the case a = in ([7]) to general a > — 1. 

2. The result in (135|) gives the Hankel transform of order |m| of the radial 
part -RIT''", and fl36l) is what one gets by inverse Hankel transformation 



of order |m|. These two results are reminiscent of, but clearly different 
from, the results in (fTOl) and ([1TD of 1481. 



3. By the asymptotics of the Bessel functions, see [T4|, 9.2.1, it is seen that 
the Fourier transform of decays as r~"~^/^ as r — )■ oo. 

5 Radon transform of generalized Zernike cir- 
cle functions; integral representation and 
DCT-formula for radial parts 

In this section, the Radon transform of the generalized circle functions 
is expressed in terms of Gegenbauer polynomials. Furthermore, an integral 
representation involving these Gegenbauer polynomials for the radial parts 
is proved, and a method of the DCT-type for computation of the radial parts 
is shown to follow from this integral representation. 

The Radon transform of is given by 

(7^Z^")(r,^)= I Z™'"(z/,/i)t// = y"i?H'"(p(t))e^-''Wrft, (46) 

Z{r,V) 
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with /(r, tp) and p, i} given in Fig. 1 for r > and < < 2tt. Thus 
p(t) = (r^+t^)^/^ , i9{t) = V^+arctan(t/r) = V^+sgn(t) arccos(r/p(t)) . (47) 
Substituting p = pit), with — oo < t < oo such that r < p(t) < 1 in the 







7 \ \ 









£(t,\|/)^\/ 





Figure 1: Line of integration 1{t,%1)) and integration point P = Q + 
t[— smifj^cosip) with — oo < t < oo, given in polar coordinates with respect 
to as (p cos p sin -(9) with p = (r^ + t^)-^/^ and ^ = ip + arctan(t/r). 



second integral in fHGl) and using the second form of d{t) in P7|) . it follows 
that 



1 

= 26^™'^ y i?|r''"(p)cos[marccos(r/p)] 



pdp 



2 _ -j-2 



= 2e^-^ / i?ri'"(p)-p^^rfp, (48) 

T 

where Tm is the Chebyshev polynomial of the first kind and of degree m, see 
[H], Ch. 22. 
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Theorem 5.1. There holds for a > — 1 and integers n, m such that 
p = I (ri — |m|) is a non- negative integer that 

(7^Z^")(r,V^) = /^±^ (l-rr^y^C:^\r)e-^ (49) 

[n + l)2a n + 2a + 1 

when < r < 1 and (7^ ZI^^''){t, tp) = for t > 1. Here C^+^ is the Gegen- 
bauer (or uhraspherical) polynomial corresponding to the weight function 
(1 - -1 < a; < 1, and of degree n, see Ch. 22. 

Proof. Consider the case that m > (the case that m < is obtained from 
this by complex conjugation). By[T3], 11.4.24, 



e-^Mt)dt = '^^;2"J-^f , c.^<l 



= , cj^ > 1 , (50) 

and so, from (HHl) 



(TIZI^'^){t,^) = e'"^^ J J e-'^'/p Jm{t)dt \ dp 



imip 'm I ^—irs 



-OO 



e-- I / K'^ip)Mps)pdp\ ds 



oo 



fmV' -n <-)a / I 1 \ f -its "^n+a+li^ ^ 



^2"(P + I)a / e—^-^^^^ds, (51) 



— oo 



where (135|) has been used, together with i™(— 1)^ = i". 

By [50], 1.12, item (10), there holds for even n {a = 1, 2n instead of n, 
u = a+l) 



oo 

o f Jn+a+l{s) , 

2 / COSTS — ds 



(52) 
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while the integral in f lS^ equals when r > 1. By [50], 2.12, item (10), there 
holds for odd n {a = 1, 2n + 1 instead of n, u = a + 1) 



-2t / sm Ts — as 



^ ^ r(n + 2a + 2) ^ ^ n v ; . 



(53) 



while the integral in f l53|) equals when r > 1. Using f l52l) and fl53|) in flSTl) . 
it follows upon some administration with and (—1)2", (— l)2{""i) that for 
< r < 1 

o2a+l 1 -p/ I 1 

{nZ-n{r,^) = {P + I)a r(n + 2a + 2) ^ ^~^'^"^'^'^"'"'^^^"^"^^ ' ^^^^ 

while {TZ ZI^'°'){t,iIj) = when r > 1. The result follows now upon some 
further administration with Pochhammer symbols. 

Theorem 5.2. There holds for a > — 1 and integer n and m such that 
n — \m\ is even and non-negative 

q + a 



)(1 -p2)-"i?H'"(p) = ^ J C^^\pcos^)e-'"'^ (I'd , (55) 



where p = \ — \m\), q = \ + \m\). 

Proof. There holds, taking |z/| = r and i/j = or n according as z/ > or 
u < in ( |49|) . with separate consideration of even m and odd m in the case 
z/ < 0, 

I Z™'"(z/,/x)d/i = K.™(l-z/2)"+i/2^«+i(^) ^ (56) 

where 

^ _ (p+1). 2^"+^r(a + l) 

'"'^"^ (n + l)2„ n + 2a + l ■ ^""'^ 

Expand 

C:^» = E /3nT'(l -p')-"^n"^'"(^,/^) (58) 
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as a function depending only on v with ^ [j? < 1. By orthogonality, see 
dSI]), there holds 



where 



^ ^ jp' + l)a 

n' + a + l {p'+ \m'\ + 1), 



Using fl56|) . it is seen that 



1 / 

-1 V-v'l^ 



^;/.'i^nw / {I - u'r^'/' c:^\u) c:;^\u) du 



where M„ follows from orthogonality of the C"'^^ as, see [16], (7.8) 

^ _ TT • 2~2«-i . r(n + 2a + 2) 
" ~ n\ (n + a + l)r2(a + l) ' 

With V = p cos -i? it then follows that 

Cri(pcos^) = 5^/3— '(l-p2)-"Zf'"(p,^) = 



irwd 



Therefore, see (126 



27r 



(59) 



(60) 



(61) 



(62) 



(63) 



^££'H2p^ - = i / Cr'(pcos^)e— '''d^ (64) 
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for integer m' such that p" := ^ — l^^'l) is a non-negative integer. An 
exphcit computation from flFTj) . f pU]) . flMl) yields now 



{p" + 1), \ -1 (/ + 1), 22°+i T{a + 1) 



n + a + 1 {p" + \m'\ + l)aJ {n + l)2a n + 2a + l 

7r-2-2""ir(n + 2a + 2) 
' n! (n + a + l)r2(a+ 1) 

r(a + i) r(g" + i)r(a + i) I g y ' ^ ' 

in which g" = p" + |m'| = ^{n + \m'\). Now replace m' by m in f lM|) and 
and results. 

Notes. 

1. Theorem 5.1 generalizes the case a = in f[T^ to general a > — 1. 
A further generalization, to orthogonal functions on spheres of general 
dimension A^ instead of disks, is provided in [51], Theorem 3.1. The 
proof of Theorem 5.1 as given here follows rather closely the approach of 
Cormack in [38]^[39] which differs from the approach used in [5T] . 

2. Theorem 5.2 generalizes the case a = of the integral representation of 
the in [ID], (A. 10) to the case of general a > — 1. Furthermore, for 
a fixed p G (0, 1), the formula f l55|) can be discretized to 

P')- Rl-'^iP) = ^ E c:^' {p -w) (66) 

k=0 

when A^ is an integer > + |m|. This yields a method of the DFT-type 
to compute the radial parts fast and reliably. 

3. The integral representation in (155|) is an excellent starting point to derive 
the asymptotics of the radial parts, by stationary phase methods, etc., 
when n,\m\ — )■ oo such that n/\m\ (0,1), with a fixed and p 
bounded away from and 1. See [S2], ENZ document. Sec. 7, item 4. 

6 Scaling theory for generalized Zernike cir- 
cle functions 
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Scaling theory for generalized Zernike circle functions is compromised by 
the occurrence of the factor (1 — p^)". In the first place, one has to restrict 
the scaling parameter e in Zl^'°'{ep,^) to the range < e < 1. While this 
restriction is quite natural, the scahng results, such as (|T7|) . for the case that 
a = allows unrestricted values of e due to polynomial form of the radial 
parts. Furthermore, the value of ^^'"(ep, -(9) at p = 1 (with < £ < 1) is in 
general finite and unequal to 0. This implies that the only natural candidate 
among the systems (Z^ '° (p, 'i9))„',m' as expansion set for Z^'"(5p, t?) is the 
case that a' = 0. Finally, an extension to shift- and- scaling theory, as in [H], 
seems also more cumbersome. Nevertheless, for the radial part there is the 
following result. 

Theorem 6.1. Let a > —1, < e < 1, and let n, m be non-negative 
integers such that p = | (n — m) is a non-negative integer. Then 

^rM= E C2ne)Rn'{p), 0<p<l, (67) 

n'=m,m+2,... 

in which the C"s can be expressed as the sum of two hypergeometric functions 
2-Fi. Furthermore, 

n' = m,m + 2, ...,n — 2 , (68) 

where p' = ^ {n — n'). Finally, when a is a non- negative integer, it holds that 

= ' n' = n + 2a + 2, n + 2a + 4, ... . (69) 

Proof. By the orthogonality condition f l3T|) for the case a = 0, there holds 
for ra' = m, m -|- 2, ... 

1 

= 2(r^' + 1) I i?r M K{p)pdp • (70) 
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Now, by ( pBj) and, subsequently (15^ with a = 0, 27rr = et, 



1 

R^'^{8p)K,{p)pdp = 



1 / oo 

\o 

oo / 1 

= {-ir2^{p+lU I I^±^^ 1 1 R^,(p)J^(ept)pdp\ dt 
\o 



(_l)P2-(p+l),(-l)("'-'")/2 [ Jn+a+lit) Jn'+l{et) ^ ^^^^ 





Using [13], first item in 9.1.27, 

^ (J„'(2;) + J„.+2(2:)) , (72) 



z 2{n' + 1) 

it is then found that 

Cn"(^) = (-l)('^+"'-'™)/'2"(p+l)„(4„,(a,£) + 4,„,+2(«,^)) , (73) 
where 

oo 

Inn"{a,£)= — dt . (74) 



Now, for n" = m,m + 2, ...,n , there holds by fl35|l 



where p" = | (ra — n"). For = n + 2, n + 4, ... , there holds by 

^-"^"'"^ = 2^r(n- + i)r(p- + c. + i) ^^'^ -"'^ +^ +''^ 

(76) 

where again p" = ^ {n — n"). When a is a non-negative integer, (175|) vanishes 
when p" + a is a negative integer, i.e., when 

n" = n + 2a + 2,n + 2a + 4, ... . (77) 



22 



When a is not an integer, infinitely many of tlie Inn" must be expected to 
be non-vanisliing. 

From f l73l) and fl75l) it is found for n' = m, m + 2, n — 2 that 



f~iin,a 



^^l)in+n'-2m)/2 2"(p + 1)^ 



-1) 



(n-n')/2 



, n — n 
2° + 1 



/ iN(n-n'-2)/2 



n — n 



(7J 



where p' = ^ (ra — n') and the proof is complete. 



7 Forward computation schemes for general- 
ized Zernike circle polynomials from ordi- 
nary ENZ and ANZ 

In this section, the generalized Zernike circle functions are expanded 
with respect to the system (Z^ )m',n' of classical circle polynomials. The az- 
imuthal dependence is in all cases through the factor exp(im'i9), and this 
allows restriction of the attention to the radial parts only. 

Theorem 7.1. Let a > —1 and let n, m be non- negative integers such 
that p = ^{n — m) is a non-negative integer. Then 



< p < 1 , 



(79) 



fc=0 



where = for A; = 0, 1, ...,p — 1 and 



Ck 



-1) 



k—p 



m + 2k + l 
m+k+p+a+1 

m + k + p + a^. _ 
m + k J 

m + 2k + l (-a)k^. 



m + p + k 
p 



a 
k — p 



(p + 1). 



m + k + p + a + 1 {k — p)\ {m + k + p + l)^ 
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^0) 



when k = p,p + 1, ... . 

Proof. By the orthogonahty condition (1311) with a = 0, there holds 



1 

a = 2{m + 2k + l) J i?™'"(p) R'^^^^Mpdp . (81) 



Using the definition fl2^ of Rn=m+2p ^m+2k^ using in the integral in 
fIST]) the substitution 

x = 2p2_i G [-1,1], p2 = |(l + x), l-p2 = |(l-x), pdp=\dx, (82) 
this becomes 



k 



1 

2{m + 2k + l) J p2™(l-p2)"pj-.™)(2p2-l)pf''")(2p2-l)pdp 



m + 2A; + 1 





1 



(1 - x)" (1 + x)'™ P|"''")(x) Pf '""^(x) rfx . (83) 



2m+a+l 

By Rodriguez' formula, see [I5], (4.3.1) or [16], p. 161, there holds 

(l-x)°(l + x)™Pp("'")(a;) = -^^^(;f)''[(l-a;)"+P(l + x)™+P] . (84) 

Then, by p partial integrations from fl83|) and ([8] 

1 

m + 2A: + l 



^ y (1 -x)^'+"(l +x)^'+'"(^)Vf™)(x)rfx . (85) 



Next, by using, see [I5], (4.2.17) or [T6], (4.14), 

^ = + « + /3 + 1) ^i-t'''^'^(^) (86) 

repeatedly, it is found that for /c > p 

dYpiO,m). , 1 (fc + m + p)! ^ 

dx) ' 2P {k + m)\ ^'-P ^""'^ 
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while this vanishes for k < p. Hence, Ck = ior k < p and for A; > p it is 
found that 

1 

m + 2k+l jk + m + py. f p{p,m-p). 



-1 



Next, again by Rodriguez' formula, 

(89) 

Thus, it is found from fISS]) and fl55]) hj k — p partial integrations that 
m + 2k + l {k + m + py {-1)''~P 



Ck 



2m+k+p+a+i [k + my p\ {k - py 

1 

-1 

m + 2k + l (A; + m + p)! (-l)^-P r(a + 1) 



2m+fe+p+Q+i [k + myp\{k-py r{a- k + p + 1) 



1 



(l-x f+P{l + x)^^"'dx . (90) 



-1 



The remaining integral in f l90|) can be evaluated in terms of F-functions as 
1 

(1 - x)°+P(l + x)^+" dx = 2™+'=+P+-+i r(fc + m + l)r(p + a + l) ^ ^^^^ 

T{k + m + p + a + l) 
-1 

The final result in (IHOj) then follows upon some further administration with 
binomials and Pochhammer symbols. 

Notes. 

1. When a is a non- negative integer, Ck vanishes for k > p + a. 

2. For the case that m = n = 0, so that p = 0, the expansion of (1 — p^)" is 
obtained, with expansion coefficients 



2fc + l fa \ , f k + a + 1 

compare [19], (10) 
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3. A short-cut of the proof of Theorem 7.1 can be obtained by noting that 
the integral on the last line of (jH3]) is essentially equal to the expansion 
coefficient Dp in the connection formula 

Pt"'\x) = j2Dr,Pir''-\x). (93) 

These connection coefficients are given in [53], Theorem 7.1.3 in terms of 
Pochhammer symbols. The proof of Theorem 7.1 as given here is "self- 
contained" , in the sense that only basic properties of Jacobi polynomials 
are used, while the proof of ^53j, Theorem 7.1.3 uses also some more 
advanced properties of the hypergeometric function 3F2. 

The result of Theorem 7.1 gives a means to transfer forward computation 
schemes from the ordinary ENZ or ANZ theory to the general setting. Below 
is an example of this. 

Theorem 7.2. Let a > —1, and let n, m be non-negative integers such 
that n — m is even and non-negative. Then the through-focus point-spread 
function U^''^{r,ip; f) corresponding to i?), see (E]), is given by 

00 

[/-'"(r,^.; /) = 27r^'"e^'"^ ^ C,V;:^,,{r,f) (94) 

A:=0 

with V^_^^2k given in semi-analytic form in f|T2|) and Ck given in f lHOj) . 
Proof. Just insert the expansion fl79p into the integral 

1 2tv 

^ifp^ ^2nipr cosi^-^) ^rn,a ^p^^)pdpdd (95) 


for f/;;^'" and use ([TT]). 
Notes. 

1. In a similar fashion, the on-axis pressure P2r((0; 0; ^) due to the ra- 
dially symmetric velocity profile v{a) = ^^^"(cr/a, i?), see f|T8l) . can be 
obtained from the on-axis pressures p2j{{0,0, z), k) due to Z2j{<j/a,0), 
given in fl20l) . and the coefficients Cj, case m = 0, in flHOl) . 

2. The availability of the through- focus point-spread functions per The- 
orem 7.2 makes it possible to do aberration retrieval, with pupil functions 
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expanded in generalized circle functions, in the same spirit as this is done 
in the ordinary ENZ theory, see Sec. |21 Similarly, radially symmetric ve- 
locity profiles, expanded into radially symmetric generalized circle func- 
tions, can be retrieved with the same approach that is used in the ordinary 
ANZ theory, see Sec. V. 

8 Acoustic quantities for baffied-piston radi- 
ation from King's integral with generalized 
circle functions as velocity profiles 

In this section various acoustic quantities that arise from baffled-piston 
radiation with a velocity profile that is expanded into generalized Zernike 
circle functions are computed in semi-analytic form. The starting point is 
King's integral, second integral expression in (fT8|) . in which V{u) is the Han- 
kel transform (|T9l) of order of v{a). Having expanded v{a) in radially 
symmetric circle functions = -R^?"' Hankel transforms 



a 









U + 1) 



a 



(aM)"+i 



J =0,1,... , 



(96) 



see fl35|) arise. 



The acoustic quantities considered are 
- pressure p((a, 0, 0) ; k) at edge of the radiator. 



- reaction force F = J pdS on the radiator, 

s 

- the power output P = f pv* dS of the radiator. 

s 

By taking z = 0, w = a in King's integral, it is seen that 



00 




(97) 







It was shown in [20], Sec. Ill from the integral result 



a 




(98) 
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by taking z = in King's integral and integrating over w, < w < a that 

oo 

F = J pdS = 2mpocka J ^^^^^^^7^ du . (99) 

S 

It was shown in [20], Sec. IV, from the representation 

oo 

p{{a,0,0),k) = ipock j ^^2^^^2)1/2 Jo{cyu)udu , < a < 00 , (100) 



as a Hankel transform and by using Parseval's formula for Hankel transforms 
that 

00 

P = j pv*dS = 2mpock j ^^^l^^.yl udu . (101) 



By inserting V = V^f into the integrals in WJ\ . fllOip . it is seen 

that the integrals 

00 







(m2 - A;2)i/2 u^+l 



du (102) 



arise. Explicitly, in ( I102p take 

- m = 0, /3 = 0, n = 2j, 7 = a for ([97D, 

- m = 0, /3 = 1, n = 2j, 7 = a for (|99|). 

- m = 2ji, (3 = a + 1, n = 2j2, 7 = « with ji < j2 for (HOip . 

These integrals will now be evaluated as a power series in ka using the method 
of [20], Appendix A. 



Theorem 8.1. Let /3 > 0, 7 > — 1 and let n, m be non-negative inte- 
gers such that n — m is non-negative and even. Then 

00 

. f Jm+fi{au) Jn+^+ijau) _ 

' J (y2_p)l/2^;3+7 - 


^ -(-l)Pa^e - (1/ + !)^ (-l/ + l-/j)p(-|/ + l),(-zfca)' 

2A;a ^ i^l + e)s T{^1 + p + + 1)T{^1 + q + 2e + 1) ' 

(103) 

where 

£ = i(/3 + 7), 5 = i(/3-7), p=i(n-m), g = i (n + m) . (104) 
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Proof. The plan of the proof is entirely the same as that of the proofs of the 
results in [20] , Appendix A, and so only the main steps with key intermediate 
results are given. 

With (m^ — /c^)^/^ as defined below (|T8i) there holds 

no 

. , Jm+^{au) Jn+^+i{au) _ 





k 







oo 

" ^ J _ ^2 



= h + h- (105) 

As to the integral Ji in f llOSp . the product of the two Bessel functions is 
written as an integral, see [54j, beginning of §13.61, and [20], (A8), from 
—ooi to ooi, the order of integration is reversed and it is used that 

h 

2g+2s+l rfo + 1 + Ffi^ 



- n2 r(g + f + s 



to obtain 

i.rri 

27ri 



1 ^^^2s+l 



(is 



Tj-s) T{2q + 2£ + 2s + 2) r(g + 1 + s)(| fcg) 

r(m + /3 + s + 1) r(n + 7 + s + 2) r(2g + 2e + s + 2) r(g + s + | 

—ooi 

(107) 

The choice of the integration contour is such that it has all poles of r(— s) on 
its right and all poles of T{2q + 2e + 2s + 2) on its left (this is possible since 
g > and /3 + 7 > —2). Closing the contour to the right, thereby enclosing 
all poles of r(— s) at s = j = 0, 1, ... with residues {—lY'^^/jl, and using the 
duplication formula of the F-function to write 

22g+2e+2i+l 

r(2g + 2e + 2j + 2) = r(g + e + j + 1) r(g + e + j + f ) , (108) 
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it follows that 

3 



j\ r(g + j + f) r(m + /?+j + l) 



, r(g + j + l) T{q + e + j + l) T{q + e + j) 
(A;a)2(9+^-)+i 



r(n + 7 + J + 2) r(2g + 2e + j + 2) 



(109) 



Replacing j + g + 1 by j = g + 1 , g + 2 , . . . , it is found after some administration 
with Pochhammer symbols (such as {x — q)q = (— 1)'^(1 — x)q) that 

J ^ -i-iya'' U + -,)e {-J + lU-J + 1 - PU-tka)'^ 

' 2ka ij+e)s r(j+p + 7 + l)r(j + g + 2£ + l) ' ^ ' 

Here it may be observed that (— j + l)q = for j = l,...,g, so that the 
summation in (11 101) could start at j = 1 as well. 

As to the integral I2 in (llOSp . again the integral representation for the 
product of two Bessel functions is used, the integration order is reversed, and 
it is used that 

J Vu^ - r(-g - s) 

k 

This yields 
/2 = i^r(i)(ia)2^ • 

T(-s) r(2g + 2£ + 2s + 2)r(-g-s-i)(Ua)2^+2'^+i 

■ as . 



T{-q - s) T{m + /3 + s + 1) r(n + 7 + s + 2) r(2g + 25 + s + 2) 

— ooi 

(112) 

The factor T{—s)/T{—q — s) is a polynomial since g = 0, 1, ... , and so the 
integrand has its poles at s = i — q — \ , j = 0, 1, ... , and at s = — r — g — e — 1, 
r = 0, 1, ... . Since e = (/3 + 7)/2 > —1, the integration contour can be chosen 
such that all poles j — g — ^ , with residues {—iy~^^/j\, lie to the right of it 
while all poles — r — g — £ — 1 lie to the left of it. Closing the contour to the 
right, and using the duplication formula again, to write 

22j+2e 

r(2j + 2e + l) = T{j + £ + i) r(j + £ + 1) , (113) 
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it follows that 



1 „• ^2e 



,r(j + l + e) r(-j + i + g) TU + e + i 



r(j + i) r(-j + i) r(j + /5 + i-p) 



(114) 



r(j+7 + f +P)r(j + 2e + | + g) 



Then some administration with Pochhammer symbols (such as r{x)/r{x — 
p) = (-l)P(l - x)p) yields 

' 2/ca ^ (j+£ + i)5 r(j+p + 7 + |)r(j + g + 2£+|) ' 

(115) 

The result in fll03p is now obtained by adding Ji in flllOp . with summation 
starting at j = 1, and I2 in flllSp . while observing that the terms j in flllOp 
yield the terms in f llOSp with even / = 2j, j = 1, 2, ... , and that the terms j 
in ffTTBj) yield the terms in ffT03|) with odd I = 2j + 1, j = 0, 1, ... . 

9 Estimating generalized velocity profiles in 
baffled-piston radiation from near-field pres- 
sure data via Weyl's formula 

In this section a brief sketch is given of how one can estimate, in the set- 
ting of baffled-piston radiation, a not necessarily radially symmetric velocity 
profile from near-field pressure data. The starting point is the Rayleigh in- 
tegral for the pressure, first integral expression in fllSp . that is written in 
normalized form as 

piu, /i ; C) = JJ v{v\ /i') ^ dv' dfj,' , (116) 

s 

where 

is the distance from the field point {u, /i, (), with C > 0, to the point {u', /i', 0) 
on the radiating surface S, for which we take the unit disk. For a fixed value 
of C > 0, the equation flll6p can be written as 

p{iy,^i;C) = {v**W{;-;0){iy,l^), (118) 
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where ** denotes 2D convolution and 
Using the Fourier transform J-", defined as 

oo oo 

{Tq){x,y)= [ [ e2™^+2"^^^g(z/,/i)di/rf/z , x,yER, (120) 



the formula (IllSp can be written as 

J'[pi;-;C)] = J'v-J'[Wi;-;C)]. (121) 

By Weyl's result on the representation of spherical waves, see [IB], Sec. 13.2.1, 
there holds 

>xn. X exp [zC((H' - (27rx)2 - (27r?/)2)i/2l 

with the same definition of the square root as the one that was used in 
connection with King's integral in f llSp . 

Now assume that the unknown velocity profile vanishes outside the unit 
disk and that it has a (1 — z^^ — yU^)"-behaviour at the edge of the unit disk, 
where a > — 1. Then v has an expansion 

t;(z/,/.) = ^Cr"^r"(^^,/^) (123) 

n,m 

in generalized Zernike circle functions, with Fourier transform 

n,m 

in which Z^^'" is given explicitly in Sec. HI Thus, when the pressure p is 
measured in the near-field plane (z/, /i ; () with ( fixed, one can estimate v on 
the level of its expansion coefficients C™'" by adopting a matching approach 
in ffmj) . using Weyl's result in ffT22D and the result of Sec. Sin (fTM]) . 
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10 Comparison with trial functions as used in 
acoustic design by Mellow and Karkkainen 



Mellow and Karkkainen are concerned with design problems in acoustic 
radiation from a resilient disk (radius a) in an infinite or finite baffle (2; = 0), 
see ^nj^lH]- The front and rear pressure distributions p+ and P-, = —p+, 
are assumed to be radially symmetric and to have the form 

00 

J2 - [cy/afr^'^ (125) 
«=o 

on the disk in accordance with the choice of trial functions used by Streng 
which is based on the work of Bouwkamp jUj. In the case that the 
normal gradient of the pressure at 2; = is considered, as is done by Mellow 
in |l5], an expansion of the form 

00 

kil - {a/afr^'^ (126) 

1=0 

on the disk has to be considered. 

In the design problem considered in [in]^[lZ], the coefficients a/ in fll25p 
are to be found such that the pressure gradient || (w, z = 0+) equals a 
desired function $(w) of the distance w of a point in the baffle plane to the 
origin. In the design problem considered in [15], the coefficients bi in (I126p 
are to be found such that p{w, z = 0+) equals a desired function \&(m;). 

The pressure p{w,z), z > can be expressed in terms of the boundary 
data p+, p_ via the dipole version of King's integral. Similarly, via the com- 
mon version of King's integral, the pressure p{w, z), z > 0, can be expressed 
in terms of the normal gradient || {w, 0) of p at z = 0. Inserting the series 
expansion (11251) and (11261) into the appropriate version of King's integral, the 
integrals 

00 

l\'±l/2 

-j Mwfi) Ji±i/2+i{afi)a dfj, (127) 


arise where a = — ^(/U^ — /c^)^/^ and where the ± follows the sign choice in 
the exponent / ± 1/2 of (1 - {(x/af) in ffT25|) and (fT26|) . To obtain ffTTfll . 
an explicit result, due to Sonine, for the Hankel transform of order of 
the functions (1 — (a/a)^)'^^/^ has been used. The integrals in (I127p are 
evaluated in the form of a double power series in ka and w/a in |45]-|47]. 
Thus, having the pressure available in this semi-analytic form, comprising 
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the coefficients a/ or bi, one can evaluate ^ {w,z = 0+) and p{w,z = 0+) 
and find the coefficients by requiring a best match with the desired function 
or '^{w). 

In the approach of the present paper, the starting point would be an 
expansion of the form 



J2ciRit'^\a/a) (128) 



1=0 



of the pressure (+-sign) or pressure gradient ( — sign) on the disk. Following 
the approach in [15]-|17], using either form of King's integral, this gives rise 
to the integrals 

oo 

//In ±1/2 
i-j Jo{wn) J2i±i/2+iian) a ^ dfi , (129) 



where now the result of Sec. H] on the Hankel transform of i?2i^^^^ has been 
used. The integral in fll29p is of the same type as the one in (11271) and can 
be evaluated by the method given in [15]-|17]. 



10.1 Numerical considerations 

In either approach, it is required to find coefficients such that a best 
match occurs between the semi-analytically computed pressure gradient or 
pressure at z = 0+, comprising the coefficients, and the desired functions $ 
or \E'. For any L = 1, 2, ... , the linear span of the function systems 

{{l-{a/afy^'/'\l = 0,l,...,L-l} (130) 

and 

{i?2f ^/^(a/a) I / = 0, 1, L - 1} (131) 

is the same. So matching using the ffist L functions in (I125p . (I126p yields the 
same result for the best matching pressure gradient or pressure at z = 0+ as 
matching using the ffist L functions in ( I128p . in theory. For small values of L, 
one finds numerically practically the same result when either system in f ll3Up . 
(I13ip is used. In the case that large values L of the number of coefficients 
to be matched are required, the approach based on (I125p . (I126p is expected 
to experience numerical problems while the one based on (I129P is likely not 
to have such problems. This is due to the fact that the functions in fll30p 
are nearly linearly dependent while the ones in fll3ip . due to orthogonality, 
are not, and this is expected to remain so after the linear transformation 
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associated with either version of King's integral Furthermore, it is to be 
expected that the semi-analytic forms, used in the matching procedure, that 
arise from any of the terms (1 — (cr/a)^)'^^'^^ must be used with much higher 
truncation levels than those that arise from the terms -R2«^^''^(cr/a). 

All this can be illustrated by comparing the expansion coefficients of 
a (1 — [a / aYY^^^"^ with respect to the system in fll3ip with those of an 
1^2^^^'^ {a / a) with respect to the system in (11301) . There is the following gen- 
eral result. 



Theorem 10.1. For m = 0, 1, ... , a > —1 and A;,p = 0, 1, ... there holds 

p-(l - p'f^- e-' = DZ:,,, ZZUp^ ^) , (132) 

1=0 

ZZ':2piP, ^) = E ^^T+2p P"^(l - pT'- , (133) 

r=0 

where 

j-^m,a _ (a + l)fc m + 2l + a + l {-k)i{m + a + l)i 

'"+2'-'= ~{m + a + l)km + k + l + a + l {a + l)/(m + k + a + l)i ' 

(134) 

j^m,a _ {a + l)p {-P)r{m + p + a + l)r . . 

(l)p (« + l).(l). • ^ ^ 

Proof. The proof of (I132p . fll34p is quite similar to the one of Theorem 7.1, 
and so only the main steps and key intermediate results are given. By or- 
thogonality, see (1ST]) , and the substitution in (15^ . there holds 

d:^:., = "'l^iLZ^' / (1 - .o'-d + pr\^) . 

(136) 

Next, Rodriguez' formula, (see flM]) ). is used, and I partial integrations are 
performed. There results for k > I 

1 



2m+k+ 



-1 

1 



4)' r(A; + /) 
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y (1 -x)^+^(l + x)'+"c/x , (137) 



-1 



while this vanishes for k < I. The remaining integral can be expressed 
in terms of F-functions as in (191]) . and then the result follows upon some 
administration with Pochhammer symbols. 

From the definition of Z^^^^p, one should find the i?'s according to 

p(.,n.)(2p2 - 1) = ^ - p^Y . (138) 

r 

By [H], 15.4.6, it holds that 

Pj"'")(2p'-1) = Pp^"''"ni-2(l-p')) = 

= ■^^^^r^F(-p,a + l + m + p; a + 1; l-p2) , 

(139) 

and then the i?'s follow from the definition of F in 114], 15.1.1. 



Example. For the case m = 0, it is found that the ratio of the r"^^ co- 

■i+2p\ 



efficient for (1 — p^Y'^'^ and the r**^ coefficient for Z^"2p(P)^) satisfy 



Dll ^ 2r + a + l (l).(a + 1). 

<2p {a + lY p + r + a + l {{p + a + lYY ' ^ ^ 

showing that the i?'s are much larger than the Z)'s. Note also that (1 — p^)^"*"" 
and Z^^2p{P) have the same L^-norm, where the weight function (1 — p^)~" 
on the unit disk is used. 



10.2 Indefinite integrals for multi-ring design 

In [47j, Mellow and Karkkainen replace the disk by a ring, and in |47] . 
Subsec. II. F, the total radiation force is expressed as an integral over this 
ring of p+ — P-- This integral can be expressed explicitly in terms of the ai 
in f ll25p since the functions (1 — (cx/a)^)^"'"^/^ have an analytic result for their 
integrals over a concentric ring. In the case that expansions involving R'^'i^^^'^ 
are used, non-trivial integrals arise. 
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Theorem 10.2. There holds for Z = 0, 1, ... 

j Rl-''\p)pdp = -^{RTip) + Rii-!2ip)) ^ (141) 

J R2l'\p)pdp = [§^^2i+2(p) + 

(4/ + 5)(4/ + l) 2' ^i") 4/ + 1 2i-2\i^)\ ' 

(142) 

where R^J:l^'^{p) = has been set for the case / = 0. 

Proof. There holds by [15j, a = in the first item in (4.1.5), 

P2i{x) = {-iyP^~'^'''\l-2x') , (143) 

where P21 denotes the Legendre polynomial of degree 21. Using this with 
a; = (1 — p'^y^'^, it is seen that 

Rir%) = (l-p^)-V^p/-/^-°)(2p^-l) = 

= (-l)'(l - p')-'/' P2z((l - p'Y^') . (144) 

For / = 0, it then follows that 

Rr'^\p)pdp = [ (i-pT'^'pdp = -(i-p^)v^ = -Rr^'ip) . (145) 



For / = 1,2,... it follows from (11441) and the substitution r = (1 — p^)^/^, 
pdp = — r dr that 



J i?°r^/'(p)prfp = (-i)'+i| PAT)dT. 

Then from [l6], (10.10), 



(146) 



Pkir) 



1 



2k + 

it follows that 



^(n+iM-n-iW) , (147) 



I Rir'^\p)pdp = (P,,+i(r) - P,Ur)) . (148) 
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Next, by [U], second item in (4.1.5), 

P2i+i{x) = {-iyxP['/'''\l-2x') . (149) 
Using this with x = (1 — p^)^/^, it is seen that 

= (-l)'P2m((l-p')'/'). (150) 

Now (141) follows from fll49p by using r = (1 - p^)^/^. 

To show (142), let / = 1, 2, ... . It follows from (11501) and the substitution 

r = (1 — p^Y^"^, pdp = —rdr that 

J R',^/\p)pdp={-ir^ j rP,,^,{r)dT . (151) 
Then by (11471) and partial integration 

j Rfl\p)pdp=^^^[TP,i^,{T)-TP,i{T)-j iP2i+2ir)-P,iir))dr . 

(152) 

Next, dg, (10.2), 

"^^'^"^^ = 2^^'-'^'^^ ^WTl^'^'^''^ ' ^^^^^ 
is used. It follows from (11521) and (11471) . (I153P and some administration that 

'-iy+^ \2l + 2 , , 41 + 3 , , 21 + 1 



41 + 3 



Al + 5 ' (4/ + 5)(4/ + l) ' Al + 1 

(154) 

Then (142) follows for / = 1, 2, ... by using ffT50|) and r = (1 -p2)i/2_ Finally, 
for Z = it follows from ffTSOj) that 

I Rr^\p)pdp = I {l-py/^pdp = -l{l-pr/\ (155) 

while (142) with / = equals by ffTMD 

- i [| P,(r) + I A(r)] = -I = -I (1 - p'f' ■ (156) 
The proof is complete. 
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